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Abstract 


Banyan  networks  comprise  a  large  class  of  networks  that  have 
been  used  for  interconnection  in  large-scale  multiprocessors  and 
telephone  switching  systems.  Regular  variants  of  Banyan  net¬ 
works,  such  as  delta  and  butterfly  networks,  have  been  used  in 
multiprocessors  such  as  the  IBM  RP3  and  the  BBN  Butterfly. 
Analysis  of  the  performance  of  Banyan  networks  has  typically 
focused  on  these  regular  variants.  We  present  a  methodology  for 
performance  analysis  of  unbuffered  Banyan  multistage  intercon¬ 
nection  networks.  The  methodology  has  two  novel  features:  it 
allows  analysis  of  networks  where  some  inputs  are  more  likely  to 
be  active  than  others,  and  allows  analysis  of  Banyan  networks  of 
arbitrary  topology.  _ 
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Introduction  j 

Banyan  networks  [2]  comprise  a  lar^e  dass  of  networks  that  have  been  used 
for  inter^nnection  in  large>scale /multiprocessors  and  telephone  switching 
systems.^  Banyan  network  is  a  network  in  which  there  is  a  unique  path 
from  each  input  to  each  output.*  Regular  variants  of  Banyan  networks,  such 
as  delta  and  butterfly  networks,  have  been  used  in  multiprocessors  such  as 
the  IBM  RP3  [6]  and  the  BBN  Butterfly  [7].  Analysis  of  the  performance  of 
Banyan  networks  has  typically  focused  on  these  regular  variants. 

Patel  [5]  presented  a  probabilistic  analysis  of  the  performance  of  delta 
networks.  His  work  assumed  that  all  sources  transmit  with  uniform  probabil¬ 
ity,  and  that  all  destinations  are  selected  with  uniform  probability.  Bhuyan 
[l]  has  extended  Patel’s  work  to  include  analysis  of  the  case  where  each  pro¬ 
cessor  has  a  single  favorite  destination  that  is  not  the  favorite  destination 
of  any  other  processor.  Kruskal  and  Snir  [3]  have  extended  Patel’s  work  by 
finding  an  asymptotic  expression  for  the  probability  that  a  destination  is 
receiving  for  networks  with  large  numbers  of  stages. 

In  what  follows,  we  present  a  methodology  for  performance  analysis  of 
general  unbuffered  Banyan  networks.  The  analysis  allows  us  to  compute  ex¬ 
actly  the  probability  of  successful  message  transmission  in  a  Banyan  network 
of  arbitrary  topology,  under  several  assumptions: 


1.  The  destination  addresses  for  messages  are  unformly  distributed  over 
the  outputs  of  the  network. 

2.  The  messages  presented  at  each  input  are  independent  of  the  mes¬ 
sages  presented  at  other  inpuliS,  and  also  of  messages  presented  on  any 
previous  cycle. 

3.  The  network  is  fully  s]mchronous,  with  all  messages  not  dropped  at 
stage  n  proceeding  simxiltaneously  to  stage  n  +  1  at  each  clock  cycle. 


Our  methodology  has  two  novel  features:  it  does  not  assume  that  all 
sources  transmit  with  equal  probability  emd  thus  allows  analysis  of  networks  or 
where  some  inputs  are  more  likely  to  be  active  than  others;  and  it  allows 
analysis  of  Banyan  networks  of  arbitrary  topology. 

Our  work  proceeds  from  the  observation  that  all  of  the  differing  topolo-  i 
gies  for  unqueued  Banyan  networks  can  be  decomposed  into  combinations  lon_ 
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of  three  basic  operatioos:  bundling,  concentration,  and  switching.  By  de¬ 
scribing  the  behavior  of  these  primitive  elements  with  a  probabilistic  model, 
we  are  able  to  evaluate  the  performance  of  any  such  network. 

We  begin  with  a  discussion  of  the  use  of  probability  mass  functions  to 
describe  network  wiring,  and  then  consider  the  effect  of  each  of  three  basic 
operations  on  these  probability  mass  functions.  Finally,  we  apply  this  model 
in  an  analysis  of  two  switching  elements,  the  common  2^  x  2^  crossbar  and 
the  Transit  RNl  switching  element. 

Modeling  Message  TrafRc  With  Probability  Mass  Functions 

A  mnlti-sts^e  network  consists  of  a  set  of  message  sources,  a  cascaded  set 
of  network  switching  elements,  and  a  set  of  message  destinations.  Often  the 
set  of  source  and  destination  nodes  is  identical. 

The  elements  comprising  a  switching  network  are  intercoimected  with 
channels.  Each  channel  consists  of  a  wire  or  group  of  wires  that  are  switched 
as  a  single  unit.  A  channel  might,  for  example,  consist  of  a  single  bidirec¬ 
tional  wire  with  serial  encoding  of  messages,  or  a  byte-wide  data  path  with 
an  associated  parity  bit. 

We  associate  with  a  channel  a  random  variable  I  whose  value  is  the  num¬ 
ber  of  messages,  or  the  load,  that  the  channel  is  carrying.  The  probability 
mass  function  (PMF)  of  this  random  variable  specifies  for  each  non-negative 
integer  j  the  probability  that  the  channel  is  carrying  j  messages.  We  call 
this  function  the  loading  probability  mass  function  (LPMF)  for  the  channel. 
If  the  random  variable  specifying  the  load  for  a  channel  is  called  I,  then  we 
denote  the  LPMF  for  the  channel  piilo). 

For  example,  a  single  channel  has  a  probability  p  of  carrying  a  message 
and  a  probability  of  1  -  p  of  being  idle.  Thus  the  LPMF  for  a  single  channel 
is  simply  the  PMF  of  a  Bernoulli  triaL 

Our  event  space  is  the  space  of  loading  configurations  for  a  particular 
network.  That  is,  if  we  define  a  network  as  a  set  of  message-carrying  wires 
connected  to  each  other  by  the  switching  elements  we  shall  define  be!  .w,  then 
the  elementary  events  in  our  event  space  are  instances  of  this  network  with 
some  load  specified  for  each  channel  in  the  network.  Obviously  the  jV  +  1 
possible  values  of  I  for  a  channel  that  can  carry  a  maximum  of  N  messages 
partition  the  event  space  into  fV  +  1  mutually  exclusive  sets  of  elementary 
events  -  each  set  containing  all  the  network  loading  comigurations  for  which 
the  load  on  that  channel  is  some  given  value. 

In  later  sections,  we  will  find  it  useful  to  associate  with  a  probability 
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mass  fimction  its  unilateral  ^-transform.  We  denote  the  ^-transform  of  a 
PMF  Pi,(*o)  by 

Bundling 

The  first  operation  is  the  simplest.  We  call  the  grouping  together  of  severzd 
fhaTin<»l«  to  form  a  single  wider  channel  bundling.  The  single  wider  channel 
that  is  a  product  of  bundling  we  sometimes  call  a  bundle.  The  loads  on  the 
constituent  rhAnnola  in  a  bimdle  must  be  independent,  as  they  will  be  in  a 
Banyan  network  with  independent  inputs.  When  we  bundle  two  channels, 
one  of  which  can  carry  between  0  and  n  messages  and  the  other  of  which  can 
carry  between  0  and  m  messages,  the  resulting  channel  can  carry  between 
0  and  n  +  m  messages.  The  loads  of  the  channels  being  bundled  are  inde¬ 
pendent  random  variables  whose  sum  we  are  forming,  so  that  the  LPMF  of 
the  resulting  bundle  will  be  the  convolution  of  the  LPMFs  of  the  component 
channels.  If  we  denote  the  bundling  of  a  and  bus  B  [pa(ao)  <  bave 

simply 


B  (p«(ao)  s  pa(ao)  •  Pfc(6o) 

where  «  denotes  convolution.  In  the  .Z-domain,  then,  bimdling  will  only 
require  forming  the  product  of  the  ^-transforms: 

Z[B\p.iao),Mbo)]]=pI{z)-d{^) 

Figure  1  depicts  the  restilt  of  bxmdling  eight  channels,  each  of  which 
carries  a  message  with  probability  1/2.  The  LPMF  is  clearly  that  of  a  bi¬ 
nomial  distribution,  because  the  sum  of  independent  identically  distributed 
Bernoulli  random  variables  is  a  binomial  random  variable. 

Concentration 

Our  second  elementary  operation  on  channels  is  called  concentration.  In 
concentration,  we  take  a  bundle  of  M  single  channels  and  form  from  it  a 
bundle  of  N  single  channels.  If  iV  <  Af,  and  the  input  bundle  is  carrying 
more  than  N  messages,  some  messages  will  be  lost. 

The  effect  on  the  LPMF  of  the  input  bundle  is  simple.  If  iV  >  Af, 
there  is  no  effect  on  the  LPMF.  If  AT  <  Af,  the  probability  that  more  than 
N  messages  can  be  carried  on  the  output  bundle  is  0,  but  in  cases  where 
messages  are  dropped,  only  enough  will  be  dropped  to  bring  the  load  to  N. 


3 


Pdoad-n} 


Figure  1:  Loading  probability  mass  function  for  an  eight-channel  bundle, 
where  each  channel  carries  a  message  with  probability  1/2. 


Thus  the  effect  of  the  operation  on  an  LPMF  pi(/o)  will  be  both  to  clip  it 
to  0  for  lo  >  N  and  to  add  to  pi(N)  the  sum  of  pt(fo)  for  /q  >  N.  Figure  2 
shows  the  result  of  concentration  on  the  LPMF  of  figure  1. 

More  explicitly,  if  the  input  LPMF  is  given  by 

<=0 

where  6{n)  is  the  unit  impulse  function,  the  result  of  iV- concentration  of 
pi(2o)i  a  bundle  composed  of  M  channels,  to  N  channels,  is  given  by 


Cm,N  [pi(^o)]  =  Pf(io)  “(iV  -  lo)  + 


j,=N+i  / 


where  u(n)  is  the  unit  step  function. 

If  the  Z-transform  of  pt(/o)  is  (2)*  then  we  have 


Z[Cm.n\pi{Io)]]  =  pfiz)  ~  ^  Pi{li)z^'  + 


f.=iV+l 


Af 

E 

ii  =  N+l 
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Figure  2:  6>concentration  of  the  LPMF  of  figure  1. 

The  first  two  terms  in  the  transform  are  the  resxilt  of  taking  the  Z~ 
transform  of  the  truncated  LPMF,  and  the  last  term  adds  in  the  ^'transform 
of  the  increased  final  element  of  the  LPMF.  Combining  the  last  two  terms, 
we  have 


2(C„j,lp,(io)ll=rf'W+  E 

li=S+l 


Switching 

The  last  operation  we  shall  be  using  is  switching.  The  switching  operation 
is  performed  on  an  input  bimdle  of  N  channels  and  specifies  the  LPMFs  for 
two  output  bundles  of  N  channels  each.  We  designate  the  output  bundles 
bundle  0  and  bundle  1.  We  specify  for  the  modeled  switch  the  probability 
(1  —  g)  that  the  each  message  in  the  input  bundle  is  switched  to  to  channel 
0;  similarly,  messages  are  switched  to  channel  1  with  probability  q. 

We  now  consider  the  LPMFs  for  the  two  output  bundles,  given  the  LPMF 
Pl(lo)  of  the  input  bundle.  Suppose  the  input  bundle  is  carrying  t  messages. 
What  is  the  probability  that  j  messages,  where  j  <  t,  will  be  switched  to 
channel  1?  It  is  the  probability  of  j  successes  in  i  Bernoulli  trials.  If  we  call 
the  random  variable  specifying  the  load  on  the  output  bundle  I,,  we  have 
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for  the  conditional  probability 


B.|l(i|i)  =  (‘V (1  - 

We  may  now  apply  the  theorem  of  total  probability  to  And  an  expression 

N 

pl.U)  = 

t=0 


where  the  summation’s  lower  botmd  is  changed  in  the  second  expression 
because  the  probability  that  the  output  channel  carries  more  messages  than 
the  input  channel  is  0. 

Thus  if  5  specifies  the  probability  of  switching  a  message  to  a  given 
output  bundle,  and  the  input  LPMF  is  given  by  pi(lo),  then  the  LPMF  for 
the  output  bundle  is  given  by 


SWIgj.Sl  2  PI.(U)  =  E  Piwf,'  V'-ll  - 


This  can  be  interpreted  as  meaning  that  the  probability  that  1,,  messages 
appear  on  an  output  bundle  is  the  probability  that  messages  were  on  the 
input  bundle  and  all  messages  were  switched  to  the  given  output  bundle, 
plus  the  sum  of  the  probability  that  /«,  -I-  1  messages  were  on  the  input 
bundle  and  exactly  of  these  were  switched  to  the  given  output  bundle, 
and  so  on,  up  to  the  maximum  possible  load  of  the  input  bundle. 

An  example  of  the  effect  of  switching  may  be  seen  in  figure  3. 

To  evaluate  the  2-transform  of  •5[pi(/o),5],  we  first  note  that  the  ran¬ 
dom  variable  describing  the  number  of  messages  on  an  output  bundle  may 
be  treated  as  the  sum  of  a  random  number  of  identically  distributed  random 
variables.  We  can  see  this  by  imagining  individually  switching  each  channel 
in  the  input  bundle  to  one  output  bundle  or  the  other,  before  considering 
whether  it  is  carrying  a  message. 

Then  there  is  one  ramdom  variable  for  each  channel  in  the  input  bundle; 
call  it  fr.  6  is  1  if  the  channel  is  switched  to  the  output  bundle  being  consid¬ 
ered,  and  0  if  the  channel  is  switched  to  the  other  output  bundle,  b’s  PMF 
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P  ( loacl»n) 

JL 

0  2  4  6  8 

Figure  3:  The  effect  of  switching  the  LPMF  of  figure  1  with  probability  0.5. 
is  then  given  by 

;%(5o)  =  (1-5)^(M  +  5^(6o-1) 

for  5  the  probability  of  switching  to  the  given  output  bundle. 

The  random  niunber  of  summands  is  the  number  of  messagi^s  that  the 
input  bimdle  was  actually  carrying;  thus  its  distribution  is  the  LPMF  of  the 
input  bundle.  To  extend  our  earlier  interpretation,  we  are  saying  here  that 
the  load  on  a  particular  output  bimdle  is  the  number  of  occupied  channels 
in  the  input  bundle  that  were  switched  to  that  output  bundle. 

Now  the  ^-transform  of  an  output  channel’s  LPMF  is  given  by  the  trans¬ 
form  of  the  sum  of  a  random  number  of  identically  distributed  random  vari¬ 
ables: 


0.3- 

0.25- 

0.2- 

0.15- 

0.1' 

0.05- 


=  pfCl  -  ^  +  Sz) 

We  note  that,  where  the  probability  of  switching  to  both  bundles  is  equal, 
the  .^-transform  for  the  LPMF  resulting  from  repeated  stages  of  switching 
has  a  particularly  simple  form.  If  b  is  the  random  variable  representing  the 
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number  of  channels  switched  to  an  output  bundle,  we  have 


^(6o)  +  6{hQ  —  1) 

2 


^  +  1 
2 


We  can  now  follow  the  rule  given  above,  so  that  n  levels  of  switching  cause 
repeated  substitutions  for  z,  and  we  have  the  recurrence  relation 


S{n) 

5(0) 


5(n-l)  +  l 
2 


with  solution 


5(n)  = 


z  +  2“  -  1 


Thus,  if  U  is  the  random  variable  for  the  load  on  the  input  bimdle  and 
In  is  the  random  variable  for  the  load  on  an  output  bundle  after  n  levels  of 
binary  switching  with  uniform  probability  of  switching  to  either  channel,  we 
have 

-1-2' 


Descriptions  of  Simple  Switching  Elements 

We  describe  two  simple  switching  elements,  the  2^  x  2^  crossbar  and  the 
Transit  RNl  switching  element,  by  using  combinations  of  our  three  prim¬ 
itive  operations:  bimdling,  concentration,  and  switching.  In  depicting  the 
primitive  operations  schematicadly,  we  use  the  symbols  shown  in  figure  4. 


The  2*  X  2*  Crossbar 

The  common  2^  x  2^  crossbar  network  is  formed  by  bundling  the  2^  inputs, 
switching  k  times  (once  per  bit  of  routing  data),  and  concentrating  the  out¬ 
puts  with  an  2^-input,  one-output  concentrator.  We  depict  the  probabilistic 
model  of  an  eight-by-eight  crossbar  in  figure  5. 

For  a  2^  X  2*  crossbar,  we  may  find  an  output  channel's  LPMF  as  follows. 
If  we  call  the  probability  that  an  input  channel  is  transmitting  Qi,  the  LPMF 
for  an  input  channel  is  given  by 

Pyivo)  =  Qi^iyo  -  1)  (1  -  <?i)^(yo) 


with  Z-trcinsform 

pj'(^)  =  +  (1  -  Qi) 
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(a)  (b)  (c) 


Figure  4:  (a)  The  symbol  for  bundling  two  input  bundles  into  one.  (b) 
The  symbol  for  concentrating  j  channels  to  k  channels,  (c)  The  s]rmbol  for 
switching  with  probability  q  to  the  top  output  channel,  and  1  —  9  to  the 
bottom  output  cbanneL 
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The  LPMF  for  a  bundle  of  2*  identical  input  channels  has  Z-transfonn 

=  (^^(■'))^ 

The  result  of  k  stages  of  switching  with  equal  probability  in  each  of  two 
directions  is 


We  note  that  this  ^-transform  is  trivially  invertible,  so  that,  setting 
Af  =  2*  and  rearranging  slightly,  we  have 


Now  we  may  perform  the  M  —*  1  concentration.  Because  this  is  a 
concentration  to  one  channel,  we  may  save  some  work  by  noting  that  we  can 
simply  consider  tl'  loading  proc.^bility  for  zero  messages  from  the  LPMF 
above;  the  concentration  forces  all  other  loading  probabilities  to  that  for  one 
message,  which  will  necessarily  be  the  complement  of  the  loading  probability 
for  zero  messages.  We  have,  for  the  loading  probability  for  zero  messi^es. 


Note  that  the  terms  in  the  summation  are  nonzero  only  where  I  =  Af,  so 
this  expression  simpUhes  to 
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P».  (0) 


The  LPMF  of  an  output  channel  is  then  given  by 

"<'«>  =  (‘  -  -  m)”) 

As  the  number  of  stages  k  in  the  crossbar  increztses,  Af  =  2*  becomes  large 
quickly,  and  pi(l)  quickly  approaches  the  limit 

-  i)  )  =  (i-'"’‘) 

In  oar  analysis,  the  probability  of  successful  message  transmission  is 
given  by  the  ratio  of  the  expected  number  of  messages  transmitted  by  all 
the  output  channels  to  the  expected  number  of  messages  on  input  channels. 
In  the  case  of  a  square  crossbar  network,  this  is  simply  pi(l)/Qi.  We  plot 
this  value  as  a  function  of  Qi,  the  input  loading  on  the  network,  for  an 
eight* by-eight  crossbar  network  in  figure  6. 

The  Transit  RNl  Switching  Element 

The  RNl  switching  element  is  a  prototype  for  the  switching  element  to  be 
used  in  the  Transit  interconnection  network  for  massively  parallel  computers, 
being  built  by  the  Transit  Group  at  MIT’s  Artificial  Intelligence  Laboratory. 
The  RNl  switching  element  can  be  configured  in  one  of  two  ways;  the  first  is 
as  two  four-by-four  crossbars,  and  the  second  is  as  an  eight-by-four  crossbar 
with  a  dilation  of  two,  meaning  that  only  four  logical  output  directions  are 
available,  but  two  messages  can  be  carried  in  each.  It  is  the  second  of  these 
configurations  whose  perform£ince  we  analyze.  We  depict  the  element  in 
figure  7. 

The  derivation  of  the  LPMF  of  a  two-channel  output  bundle  for  the  RNl 
switching  element  follows.  As  above,  if  we  call  the  probability  that  an  input 
channel  is  transmitting  Qi,  the  LPMF  for  an  input  channel  is  given  by 

Py(yo)  =  QiHvo  -  1)  +  (1  -  C?i)^(yo) 


ll 


P( Success) 


Figure  6:  The  probability  of  successful  message  transmission  as  a  function  of 
Qi,  the  input  loading  on  the  network,  for  an  eight-by-eight  crossbar  network. 


8->2 


Figure  7:  The  RNl  switching  element,  in  the  eight- by- four,  dilation  two 
configuration. 
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with  2^- transform 


Pv(^)  =  <3*-'  +  (l-<?i) 

The  LPMF  for  a  bundle  of  8  identical  input  channels  has  ^-transform 

The  resiilt  of  two  stages  of  switching  with  equal  probability  in  each  of  two 
directions  is 

=  («.  (i^).  a -«,))• 

=  i:e)(-w(^r 

Again  as  in  the  case  of  the  crossbar,  we  invert  the  transform 

-•<*-)  =  (f)’(|:0(^ 

and  then  perform  the  concentration.  In  this  case  the  concentration  is  to 
two  channels,  so  that  we  must  consider  probabilities  for  the  two  cases  that 
the  output  bundle  carries  zero  messages  or  one  message  in  order  to  use  the 
method  we  did  before  for  deriving  the  concentrated  LPMF. 

For  zero  messages,  the  sum  is  zero  except  where  /  =  8,  so  that  we  have: 

For  one  message,  the  stim  is  zero  except  where  1  =  7,  so  we  have: 


=  20, 
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After  concentration,  the  probability  for  two  messages  must  be  the  sum 
of  the  probabilities  for  higher  loads,  so  that  the  LPMF  for  a  two-channel 
output  bundle  is  given  by 


p,(/o)=  (l-^)%(io)  + 


+ 


S(lo  -  2) 


We  form  the  probability  of  successful  message  transmission  as  the  ratio 
of  the  expectation  of  the  number  of  messages  on  all  output  channels  to 
the  expectation  of  the  number  of  messages  on  all  input  channels.  In  this 
analysis  we  have  assumed  uniformity  and  independence  of  input  loading  and 
a  uniform  distribution  of  message  destinations,  so  that  the  expectation  of 
the  input  loading  is  simply  12n=i  Qi  =  ^Qi  output  loading 

(if  we  recall  that  the  random  variable  giving  the  number  of  messages  on  an 
output  bundle  is  1)  is 


EM  =  ■*  (i  ■  (i  -  t)V  '•('-('-  t)  '('  +  ^))  ) 


Thus  the  probability  of  successful  message  transmission  is  given  by 


PsMT 


Qi 

Qi 

1_(i  +  32,) 

Qi 


We  plot  the  probability  of  successful  message  transmission  versus  the 
input  loading  in  figure  8. 


Analyzing  the  Performance  of  More  Complex  Networks 

It  may  be  difficult  to  simplify  the  expressions  describing  more  complex  net¬ 
works  built  from  arrays  of  simple  switching  elements  like  those  we  have 
analyzed  above.  Indeed,  Patel  [5)  and  Kruskal  and  Snir  [3]  derive  expres¬ 
sions  only  for  simple,  regular  networks;  these  are  o'*  x  6**  delta  networks  in 
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Figure  8:  Probability  of  successful  message  transmission  plotted  vs.  input 
loading  for  the  Transit  RNl  switching  element  in  its  eight- by-four,  dilation 
two  configuration. 

the  case  of  Patel’s  work,  and  square  Banyan  networks  of  arbitrary  dilation 
in  the  case  of  Kruskal  and  Snir’s  work. 

The  advantage  of  our  approach  is  that  such  analyses  are  automated.  One 
specifies  the  topology  of  the  network,  forms  the  sequence  of  operators  that 
describes  the  loading  probability  mass  function  for  an  output  bundle,  and 
evaluates  it.  This  evaluation  can  be  numeric  or  in  the  form  of  a  parameter¬ 
ized  expression.  The  derivation  of  such  an  expression  for  a  complex  network 
is  aided  by  the  use  of  a  symbolic  mathematics  program  like  Maesyma  or 
Mathematica.  We  present  a  set  of  Mathematica  functions  that  may  be  used 
for  such  analysis  in  the  appendix. 

Future  Work 

The  methodology  described  above,  despite  its  advantages,  does  not  yield 
a  completely  satisfactory  model  of  a  Bsmyan  network’s  performance.  We 
describe  now  some  of  the  disadvantages  of  the  methodology. 

The  probability  of  successful  message  transmission  alone  will  not  be  a 
faithful  meaisure  of  performance  in  a  network  that  is  buffered.  In  the  Transit 
network,  for  example,  eilthough  the  individual  switching  elements  themselves 
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do  not  contain  buffers,  messages  are  effectively  buffered  at  the  inputs  to  the 
network.  Thus  it  will  be  desirable  to  extend  the  methodology  with  queueing- 
theoretic  techniques  to  allow  the  creation  of  more  faithful  models.  In  future 
work,  we  hope  to  do  this  in  a  fashion  that  continues  to  allow  analyses  of 
Banyan  networks  of  arbitrary  topology. 

While  the  methodology  does  allow  analysis  of  networks  where  one  or 
more  sources  are  likely  to  be  more  active  than  others,  it  does  not  easily  lend 
itself  to  an  analysis  of  networks  where  one  or  more  sinks  are  more  likely  to 
be  the  destination  of  messages  than  others.  In  fact,  the  general  case  of  this 
problem,  where  messages  entering  a  Banyan  network  of  arbitrary  topology 
have  an  arbitrary  distribution  of  destination  addresses,  remains  imsolved  to 
date. 

The  analysis  technique  described  is  appropriate  only  to  Banyan  net¬ 
works.  While  these  constitute  a  large  class,  some  of  the  fault-tolerance 
features  of  Banyan  networks  used  in  practice  may  include  redundant  paths 
between  sources  and  sinks.  It  will  be  necessary  to  extend  the  technique  to 
encompass  replications  and  dilations  of  Banyan  networks;  in  more  compli¬ 
cated  cases,  it  may  be  necessary  to  supplement  it  with  a  different  approach, 
or  abandon  it  altogether. 

Another  disadvantage  of  the  methodology  we  have  described  lies  in  its 
tacit  assumption  that  the  network  modeled  is  completely  synchronous.  This 
assumption  is  not  always  justified;  for  example,  in  the  case  of  a  circuit- 
switched  network  like  the  Transit  network,  successful  message  transmission 
creates  a  circuit  which  is  held  open  until  a  reply  is  sent.  The  circuit  is  held  for 
a  number  of  cycles,  during  which  other  messages  may  be  transmitted  from 
the  inputs  and  be  blocked  because  paths  at  succeeding  stages  are  already  in 


use. 

A  related  disadvantage  of  our  methodology  lies  in  its  asstimption  that 
the  path  being  built  by  a  message  being  transmitted  in  a  circuit-switched 
network  immediately  disappears,  freeing  all  associated  resources,  if  the  mes¬ 
sage  is  blocked,  whereas  in  reality  it  will  take  a  number  of  cycles  for  these 
resources  to  be  freed.  Nussbaum  and  his  colleagues  have  found  this  to  be  a 
significant  factor  in  discrepancies  between  Patel’s  model  and  their  simula¬ 
tion,  as  described  in  in  [4]. 

We  hope  to  address  some  of  these  disadvantages  in  ongoing  work.  The 
ideal  result  would  be  a  technique  for  automatically  generating  an  accurate 
model  of  the  performance  of  a  multistage  interconnection  network  given  only 
a  description  of  the  network  topology. 
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Appendix:  Mathematica  Functions  for  Banyan  Network  Anal¬ 
ysis 

concentrate: : usage  a 

"concentrate Cx,  n]  concentrates  the  LPHF  x  to  n  channels." 
concentrate Cx_,  n_] 

(*  get  distribution  lor  0  through  u-1  chauuiels ,  and  add 
as  last  element  the  sum  ol  the  rest  of  the  channels.  «} 
Append [Take Cx.  n] .  Apply [Plus.  DropCx,  n]]] 

discreteconvolution: : usage  > 

"discreteconTolutionCx.  y]  treats  x  and  y  as  0-based 
vectors  and  returns  their  discrete  convolution." 

discreteconvolutionCx..  y_]  :* 

Block  C-Cxlgth,  ylgth,  Igth}. 
xlgth  •  Length Cx] ; 
ylgth  ■  LengthCy]; 

Igth  ■  xlgth  ylgth  -  1; 

(«  in  summation,  portions  ol  sequence  with  Indices 
out  ol  range  lor  sequences  must  be  treated  as 
0.  *) 

Table CSumCllCk  <1  II  k  >  xlgth  II 

(n-k-^1)  <111  (n-k-^1)  >  ylgth, 

0. 

(*  because  ol  the  0->l  index 

translation,  se  increase  the  y-index 
to  shilt  the  result  sequence  back 
doun  to  begin  at  1.  «) 
xCCk]]  yCCn-k+1]]], 

•Ck ,  xlgth}] , 

{n,  Igth}]] 

bundle : : usage  * 

"bundle Cx,  y]  lorms  the  LPNF  that  results  Irom  bundling 
tvo  input  bundles  sith  LPNPs  x  and  y." 

bundle Cx_,  y_]  :« 

discreteconvolutionCx,  y] 
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switch: : usage  > 

"switch Cz,  p]  ratuzns  the  LPMF  of  an  output  bundle  to 
which  z  is  switched  with  probability  p." 

switch Cz..  p_]  := 

Block [{Igth}, 

Igth  «  Length Cz] ; 

TableCSufflCzCCi'^1]]  BinomialCi,  n]  p*n  (l-p)“(i-n), 
■Ci,  n,  lgth-1}], 

{n,  0.  lgth-1}]] 
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